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ABSTRACT 

Context. In the past decade, more than one hundred asteroid models were derived using the lightcurve inversion method. Measured 
by the number of derived models, lightcurve inversion has become the leading method for asteroid shape determination. 
Aims. Tens of thousands of sparse-in-time lightcurves from astrometric projects are publicly available. We investigate these data 
and use them in the lightcurve inversion method to derive new asteroid models. By having a greater number of models with known 
physical properties, we can gain a better insight into the nature of individual objects and into the whole asteroid population. 
Methods. We use sparse photometry from selected observatories from the AstDyS database (Asteroids - Dynamic Site), either alone 
or in combination with dense lightcurves, to determine new asteroid models by the lightcurve inversion method. We investigate various 
correlations between several asteroid parameters and characteristics such as the rotational state and diameter or family membership. 
We focus on the distribution of ecliptic latitudes of pole directions. We create a synthetic uniform distribution of latitudes, compute 
the method bias, and compare the results with the distribution of known models. We also construct a model for the long-term evolution 
of spins. 

Results. We present 80 new asteroid models derived from combined data sets where sparse photometry is taken from the AstDyS 
database and dense lightcurves are from the Uppsala Asteroid Photometric Catalogue (UAPC) and from several individual observers. 
For 18 asteroids, we present updated shape solutions based on new photometric data. For another 30 asteroids we present their partial 
models, i.e., an accurate period value and an estimate of the ecliptic latitude of the pole. The addition of new models increases the 
total number of models derived by the lightcurve inversion method to ~200. We also present a simple statistical analysis of physical 
properties of asteroids where we look for possible correlations between various physical parameters with an emphasis on the spin 
vector. We present the observed and de-biased distributions of ecliptic latitudes with respect to different size ranges of asteroids as 
well as a simple theoretical model of the latitude distribution and then compare its predictions with the observed distributions. From 
this analysis we find that the latitude distribution of small asteroids (D < 30 km) is clustered towards ecliptic poles and can be 
explained by the YORP thermal effect while the latitude distribution of larger asteroids (D > 60 km) exhibits an evident excess of 
prograde rotators, probably of primordial origin. 

Key words, minor planets, asteroids - photometry - models 



1. Introduction 

The lightcurve inversion method (LI) is a powerful tool that al- 
lows us to derive basic physical properties of asteroids (the ro- 
tational s tate and the shape) from their disk-integrated photom- 
etry ( see iKaasalainen & Torppall200Tl ; iKaasalainen et al.ll2001l 
I2002I) . This photometry can be dense-in-time, which typically 
consists of tens to a few hundreds of individual data points ob- 
served during one revolution. This is in contrast to sparse-in- 
time, where the typical separation of individual measurements is 



large compared to the rotation period. For sparse data, we usu- 
ally have a few measurements per night, such as in the case of 
astrometric sky surveys. In the following text, we use the terms 
"dense lightcurves" and "sparse lightcurves". 

To obtain a unique spin and shape solution, we need a set of 
at least a few tens of dense lightcurves observed during at least 
three apparitions. Based on simulated photometric data and the 
survey cadence of the Panoramic Survey Telescope an d Rapid 
Response System (Pan-STARRS), IKaasalainen] (2004) showed 
that we can also use only sparse data for the inversion technique. 
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In this case, a unique model can be derived from more than about 
one hundred calibrated measurements observed during 3-5 years 
if the photo metric accuracy is better than ~5% dDurech et al.l 
I2005L 120071) . Sparse data available so far are not that accurate. 
Nevertheless, for many asteroids with high lightcurve ampli- 
tudes, it is possible to derive their models from current sparse 
data (see Section[2]for more details). We can also combine sparse 
and dense data to derive models. First results from this approach 
were presented bv lDurech et al.1 (120091) . where sparse data from 
the US Naval Observatory in Flagstaff (USNO) were used. 

Currently (January 2011), there are 113 models of aster- 
oids derived by the lightcurve inversion method; most of them 
are stored in the Data base of Asteroid M odels from Inversion 
Techniques (DAMIT 1 . 15urech et alJfeOlOl) . Most of these mod- 
els were derived from dense lightcurves. Only 24 of them were 
comp uted from combined dense and sparse data dDurech et al.l 
120091) . The AstDyS database (Asteroids - Dynamic Site 2 ), 
which contains data from astrometric projects, is another pos- 
sible source of sparse data. However, most of the data are 
not accurate enough to be used for inversion alone. On the 
other hand, even noisy sparse data in combination with a few 
dense lightcurves can give us, in many cases, a unique solution 
dDurech et al.ll2007h . The aim of our work was to gather these 
data, keep only those that were useful, and then combine them 
with dense lightcurves in the lightcurve inversion method. 

Dense data are best used to define the rotational period 
and constrain the period interval that must be searched during 
the model computation (see Section 13.11 for more details). On 
the other hand, sparse data usually cover a long time interval, 
typically over several apparitions, and carry information about 
brightness variations for different geometries, which constrains 
the pole directions. 

A priori information about rotational periods of asteroids 
plays an important role in the process of model determination. 
When an approximate period is known, we search for the solu- 
tion near this value (details in Section I3.lt and thus save con- 
siderable computational time. We use the latest update of the 
Minor Planet Lightcurve Database 3 published by I Warner et al.l 
d2009l) to check for previously derived periods. For many aster- 
oids, there are only a few sparse lightcurves from different astro- 
metric observatories available but no dense lightcurves. In these 
cases, we must scan the whole interval of expected period values 
(2-30 hours). This approach is time-consuming and there is no 
guarantee that the correct period will lie in the scanned interval. 

The knowledge of rotational states of asteroids is fundamen- 
tal for understanding the history of the Solar System, specifically 
the accretion of planets or the collisional processes. For exam- 
ple, it was presumed that due to collisional evolution, the spin- 
vector distribution of Main Belt Asteroids (MBAs) should be 
nearly isotropic, possibly with a small excess o f prograde spins 
dDavis et alJll989h . lJohansen & Lacerdald2010l) performed a hy- 
drodynamical simulation of the accretion of pebbles and rocks 
onto protoplanets and speculated that the trend of prograde rota- 
tors among the largest asteroids is primordial. 

First statistical analys es of the s pin-vec tor distribu- 

tion were presented by | Magnussqn| (1 19861 1 19901) and 
Drummond et al.l (1 19881 1 19911) . later bv lPravec et al.l (|2002|) and 



Skoglov & EriksorJ d2002[). They all observed a lac k of poles 
~ el - 



close to the ecliptic plane. iKrvszczvriska et al.l d2007l) used more 



objects in the analysis, finding that the distribution was strongly 
anisotropic with a moderate excess of prograde spins in the lim- 
ited size range from 100 to 150 km. Interpretation of this de- 
population of poles close to the ecliptic plane is still unclear. 
Probable candidate s are selection effects, the role of inclination , 
the YORP effect 4 (iRubincaml 120001 1 Vota-ouhlickfetai] 120031 
or a combination of these. The YORP effect acts only on small 
bodies with D < 40 km. Asteroi ds with these sizes hav e non- 
Maxwellian spin rate distribution dPravec & Harris || 2000l) and is 
partic ularly evident for asteroids with D < 14 km d Warner et al.l 
120091) . It is believed that the YORP effect is responsible for this 
trend since it can either spin up or spin down an irregularly- 
shaped asteroid on the timescale shorter than the typical time 
between collisions and also affects the obliquity of spin axes 
dRubincamll2"000l:lBottke eTal]|2006l) . 

In the Near Earth Asteroids (NEAs) population, the lat- 
itude distributi o n of poles is different from th at of MBAs 
dLa Spina et al.l l2004t IKrvszczvriska et al.1 120071) . i.e., there 
is a significant excess of retrograde spins probably caused 
by the transport mechanism of MBAs to Earth-crossing 
space by gravitational resonance s and the Yarkovsky effect 5 
dMorbidelli & Vokrouhhckv1l2003l) . There is no statistically sig- 
nificant ck^tering_jnjhetongjtude of poles of either MBAs or 
NEAs dKrvszczvriska et al.ll2007l) . 

As the number of asteroid models with known physical prop- 
erties grows, we can have a better insight into the nature of in- 
dividual objects and into the asteroid population as a whole. In 
Section [2] we describe available dense and sparse photometric 
data and the selection of astrometric observatories with qual- 
ity sparse data. In Section [3] we present new asteroid models 
derived from combined photometric data sets or from sparse 
data alone, mentioning a few individual objects and define sev- 
eral procedures on how to test the reliability of new mod- 
els. In Section |U we present a statistical analysis of asteroid 
physical parameters that we derived using the lightcurve in- 
version method or adopted from different sources (proper ele- 
ments from the AstDyS database, diameters from IRAS, . . . ). 
We also present results of a numerical simulation that allowed 
us to estimate the bias in pole directions of the lightcurve inver- 
sion method. Using these results, we then corrected the observed 
pole distributions for this effect. Finally, in order to explain the 
observed latitude distributions, we present a simple theoretical 
model of the latitude distribution in Section[5] 



2. Photometric data 

The main source of dense photometric lig htcurves is the Uppsala 
Asteroid Photometri c Catalogue (UAPC, lLagerkvist et al.ll987t : 
iPiironen et al 1 120011) . where the lightcurves for about 1 000 as- 
teroids are stored. We also used data from several individual ob- 
servers (Table [6). 

Sparse photometry was first used in combination with dense 
data for lightcurve inversion by IDurech et al.l d2009l) . These 
sparse data were from the USNO-Flagstaff station and had a 
typical photometric uncertainty of ~8— 10%. Other sparse pho- 
tometric measurements are produced by many astrometric sur- 
veys, but mostly as a by-product. In most cases, asteroid magni- 
tudes are given to only one decimal place, i.e., the accuracy is 0. 1 
mag at best. Whether or nor this is sufficient for a unique shape 



1 http : //astro . troja.mff . cuni . cz/projects/asteroids3D 

2 http://hamilton.dm.unipi.it/ 

3 http : // cf a-www . harvard . edu/ iau/lists/Lightcurve- 
Dat.html 



4 Yarkovsky-O'Keefe-Radzievskii-Paddack effect, a torque caused 
by the recoil force from anisotropic thermal emission 

5 a thermal force acting on a rotating asteroid 
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determination for reasonable number of asteroids can be de- 
duced from asteroids lightcurve amplitude distribution. We used 
lightcurve amplitude data fo r ~2 500 asteroids from the Minor 
Planet Lightcurve Database dWarner et al.ll200 9) and found that 
the mean lightcurve amplitude is ~0.3 mag. For 19% of aster- 
oids, the amplitude is >0.5 mag. This means that, in principle, 
photometry with an accuracy of ~0. 1 mag carries sufficient in- 
formation about rotational states and shapes for a significant 
number of asteroids. 

Our goal was to find out which observatories produce pho- 
tometry suitable for lightcurve inversion and to use these data for 
determining new asteroid models. Through to September 2009 
(the time of the data download), data for more than 350000 
objects from almost 1 500 observatories were archived on the 
AstDyS server. Some of the observatories contributed with only 
a few data points, while others contributed tens of thousands of 
photometric measurements (e.g., large sky surveys such as the 
Catalina Sky Survey, LONEOS, or Siding Spring Survey). 

2.1. Data reduction 

The quality of the sparse photometry archived on the AstDyS 
varies significantly. We investigated the photometry carefully by 
establishing criteria for its quality. Then, using those criteria, we 
choose only those data that were useful for inversion. 

For each observatory, we extracted photometric data for the 
first 10000 numbered asteroids if there were at least 30 data 
points for a single lightcurve. We then transformed this pho- 
tometry to the standard format used in lightcurve inversion: we 
computed geometry of observation (astrocentric ecliptic coordi- 
nates of the Sun and the Earth), corrected for light-time, normal- 
ized the brightness to the distance of 1 AU from the Sun and the 
Earth, and excluded clear outliers. 

For further investigation, we selected 13 observatories that 
fulfilled the condition of having data for more than ~50 as- 
teroids. This resulted in almost 30000 sparse lightcurves for 
~9 000 asteroids. In the next step, we estimated mean uncertain- 
ties of individual observatories and, based on these uncertainties, 
we assigned a relative weight to the data from each observatory. 
In this process, we assumed that the brightness vs. solar phase 
angle relation can be fitted with a simple relation for each sparse 
lightcurve: 



f(a) = cos 2 ^— j a exp (— — j + c a + d 



(1) 



where a is the solar phase angle 6 and a, b, c and d are free pa- 
rameters. Then, we constructed a histogram of residuals (rms) 
for each observatory comparing actual data against the model 
given by Eq. Q] Four examples are plotted in Fig. Q] The dis- 
persion is caused by observational uncertainties and by the am- 
plitudes of the lightcurves. From these histograms, we esti- 
mated the "FWHM" 7 values and the most frequent residual (the 
mode); median values of the residual distributions for each lo- 
cation (Table [1]) were computed. Observatories with a high me- 
dian or "FWHM" value (>0.2 mag for both) are not suitable for 
the lightcurve inversion (e.g. observatory 691 in Fig. [1}. Data 
from only seven observatories, listed in Table Q] with non-zero 
weights, had sufficient accuracy and so could be used for mod- 
eling. Based on the values of medians and "FWHMs", we esti- 
mated a weight for the photometric data from each observatory 



500 
400 
300 
200 



0.1 0.2 0.3 0.4 0.5 
rms 




6 the Sun-asteroid-Earth angle 

7 the width of the distribution in the half of its maximum 



Fig. 1. Four histograms of residuals comparing actual data 
against the model given by EqQ] of all sparse lightcurves be- 
longing to the following observatories: 689 USNO, 691 Steward 
Observatory, 703 Catalina Sky Survey and E12 Siding Spring 
Survey. Number of bins is ~ y/N, where N is the total number 
of sparse lightcurves used for histogram construction. 



relative to dense data, which has a unity weight (see Table[TJ- We 
assumed that the typical accuracy of dense lightcurves is ~0.02 
mag. 

The USNO in Flagstaff (MPC code 689) and the Hipparcos 
satellite are clearly the best observatories with respect to photo- 
metric accuracy. Other observatories are less accurate but they 
still hold some information about rotational states and shapes of 
asteroids. For any given asteroid, we have typically 2^1 sparse 
lightcurves from different observatories covering the last ~10- 
15 yea rs. Data from USNO were already used by IDurech et al.l 
(2009). We updated those data along with adding another 20- 
30% of new data from the years 2008 and 2009 if there was an 
apparition for the asteroid. 



3. Results 

3. 1 . Models of asteroids 

Our lightcurve inversion (LI) method is based on the opti- 
mization of unknown parameters of the shape (modeled as a 
convex hull), the rotational state, and the scattering law (see 
iKaasalainen & Torppall200ll iKaasalainen et alJl200lh . The pa- 
rameter space have many local minims. Since LI is a gradient- 
based method that converges to a local minimum near the ini- 
tial choice of parameter values, it is critical to find the global 
minimum for the data set and then do the modeling. Finding a 
global minimum involves a systematic search through all rele- 
vant parameter values. Each model corresponding to a particular 
local minimum is characterized by a single value x 2 , which cor- 
responds to the quality of the fit. 

A unique solution is defined as follows: (i) the best period 
has at least 10% lower x 2 than all other periods in the scanned 
interval, (ii) for this period, there is only one pole solution with 
at least 10% lower x 2 than the others (with a possible ambiguity 
for A + 180°), and (iii) this solution fulfills our additional tests 
(see Section l3~3l . 
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Table 1. Comparison of estimated characteristics of residuals for 13 selected observatories: mode, "FWHM" and median. 



Obs 


A'lc 


Mode 


FWHM 


Median 


N P 


Weight 


Observatory name 


608 


2459 


0.26 


0.20 


0.27 


37 





Haleakala-AMOS 


644 


2567 


0.22 


0.22 


0.24 


36 





Palomar Mountain/NEAT 


683 


218 


0.18 


0.25 


0.20 


39 





Goodricke-Pigott Observatory, Tucson 


689 


1970 


0.14 


0.12 


0.15 


118 


0.3 


U.S. Naval Observatory, Flagstaff 


691 


1893 


0.23 


0.22 


0.24 


39 





Steward Observatory, Kitt Peak-Spacewatch 


699 


546 


0.17 


0.11 


0.18 


33 


0.1 


Lowell Observatory-LONEOS 


703 


8350 


0.17 


0.16 


0.19 


54 


0.15 


Catalina Sky Survey 


704 


8333 


0.42 


0.17 


0.42 


311 





Lincoln Laboratory ETS, New Mexico 


950 


80 


0.14 


0.11 


0.15 


180 


0.15 


La Palma 


E12 


1354 


0.14 


0.16 


0.15 


41 


0.1 


Siding Spring Survey 


G96 


1810 


0.14 


0.22 


0.17 


43 


0.1 


Mt. Lemmon Survey 


H07 


161 


0.20 


0.18 


0.23 


47 





7300 Observatory, Cloudcroft 


Hip 


49 


0.10 


0.10 


0.11 


53 


0.3 


Hipparcos satellite 



Notes. For each observatory, the table gives also the number of sparse lightcurves A'lc, and the average number of data points for a single 
lightcurve Np. 



The most time-consuming part of the lightcurve inversion 
method is scanning through all periods within a chosen interval, 
which we determined by using the periods and reliability codes 
given in the Minor Planet Lightcurve Database 8 . For each aster- 
oid, we scanned an interval centered at the reported period value 
P with a range of + 1%, +5% and +20% of P for reliability codes 
4, 3 and 2, respectively. Half and double periods were tested later 
on. 

We combined relative lightcurves from the UAPC and from 
individual observers together with sparse data obtained from the 
AstDyS site to create a data set for each asteroid. This gave us 
data sets for ~2 300 asteroids (in ~900 cases there were only 
sparse data available) to which we then applied the lightcurve 
inversion method and then we ran the additional tests described 
in Section 1331 We derived 80 new unique models, 16 of which 
are based only on sparse data. Basic characteristics of these mod- 
els are listed in Table[3] We estimated the uncertainty in the pole 
direction as ±10-20° based on previous results with limited data 
sets. As might be expected, the uncertainty seems dependent on 
the number of dense and sparse photometric data. The longitude 
uncertainty increases for higher latitudes because meridians on 
a (A, /3)-sphere are more dense with increasing latitude. These 
uncertainties are discussed in more detail in Section 14.31 The 
uncertainty of the rotational period depends on the time inter- 
val covered by observational data and is of the order of the last 
decimal place of period values P in Tables [3] [4] and [5] 

In some cases, we were able to determine a unique rota- 
tional period, but we had multiple pole solutions with similar 
ecliptic latitudes /3. These models give us accurate period val- 
ues and rough estimates of ecliptic latitudes j3, which are also 
important parameters. In Table |4] we present results for 30 par- 
tial models, where (3 is the mean value for all different models 
if the dispersion is smaller than 50°. We defined a parameter 
A = |/?max-/8minl/2 as being the estimated uncertainty of /3, where 
Armx and Anin are the extremal values within all /3. 

All new unique shape models are now included in DAMIT. 



3.2. Comments to selected models 

In DAMIT there are several solutions designated as "prelimi- 
nary". These models do not have a well-constrained pole so- 
lution or are based on combined data sets. For 18 of those as- 
teroids we derived updated model solutions based on additional 
photometric data (see Table 0. The difference between the old 
and new model for asteroid (1223) Neckar was significant. The 
new model has a slightly different period, but the pole direc- 
tions and shapes are nearly similar to the old model. The current 
data suggest a period of P = 7.82401 hours (previous value was 
P = 7.82123 hours). 

The asteroid (4483) Petofi was recently observed by Brian 
Warner. We derived a shape solution from three poor dense 
lightcurves and one sparse lightcurve from Catalina Sky Survey. 
Warner used these four lightcurves in combination with his new 
observ ations and also derived the shape model of Petofi dWarnerl 
1201 lah . His period of P = 4.3330 hours and pole direction (90°, 
35°) are close to our solution of P = 4.33299 hours and (107°, 
40°). 

The asteroi d (83 2) Karin was also studied by 
ISlivan & Molnarl (1201 Oh : their solution with P = 18.352 
hours and pole (51° or 228°, 41°) confirms our results, see Table 

m 

In past decades, occultations of stars by several asteroids 
were observed. These events give us additional information 
about the shape and can help resolve which mirror solutions 
of a model is the correct one. According to the recent work of 
IDurech et al.1 (1201 ll) . asteroid occultation measurements prefer 
pole solutions of (122°, -44°) for (10) Hygiea, (347°, 47°) for 
(152) Atala, (28°, -72°) for (302) Clarissa, (223°, 67°) for (471) 
Papagena, and (296°, 41°) for (925) Alphonsina. Spin solutions 
preferred by asteroid occultations appear in bold font in Tables 
E]and[5] 



3.3. Models and method testing 



We constructed five additional tests to be sure that the new mod- 

s |http : //cf a-www . harvard . edu/iau/lists/LightcurveDat . htSfl l are reliable. We performed the first two tests for all models, 
there is also an explanation and more details about the reliability codes For models derived only from sparse data, which are presented 
too for the first time, we performed three additional tests: 
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Inertia tensor of the shape model. The lightcurve inversion 
method we use assumes that asteroids are in a relaxed rotational 
state, which means that derived models should rotate around the 
axis with a maximum moment of inertia. For each derived shape, 
we computed princi pal moments using equations presented by 
iDobrovolskisld 19961) and checked if the rotation axis was close 
to the principal axis of the maximum momentum of inertia. We 
rejected those models for which the angle between the spin axis 
and the axis with a maximum momentum of inertia was larger 
than +30°. However, this criterion is too strict for elongated 
models with similar sizes along the rotational axis and the axis 
that is both perpendicular to the rotational axis and is the min- 
imal size of the model. In this case, the principal moments for 
these two axes are similar. Under these circumstances, the angle 
between the spin axis and the axis with a maximum momentum 
of inertia can be large even for realistic shapes and so we allowed 
these models to pass this test. 

Half- and double-period models. In cases where we have only 
a few dense lightcurves for a given asteroid, it is easy to con- 
fuse the correct rotational period with its half or double value. 
When the a priori period value in the Minor Planet Lightcurve 
Database is uncertain, which corresponds to a low reliability 
code (see Section 13. U . it is reasonable to check if the half and 
double period value give a better fit. If the period was in doubt, 
we searched for a solution also around 2P and P/2; if the^ 2 was 
lower than I Ax 2 of the solution with period P, we rejected the 
model as unreliable. 

Reduction of the number of sparse photometric data. In this 
test, we used only the sparse data sets for modeling. For 63 as- 
teroids, this led to unique shape and spin state solutions (after 
performing the tests described above). For each solution, we 
randomly reduced the original amount of observed sparse data 
points to 90% and used these new limited data sets again in the 
lightcurve inversion. Our expectation was that we would not get 
a unique solution when using less data. This was true for five as- 
teroids. These models, when using the full sparse data sets, are 
not necessarily wrong, but the amount of available data is prob- 
ably just at the level when a unique solution can be derived. The 
important point of this test is that, for a given asteroid, we did not 
find two different but formally correct solutions when using the 
full versus reduced data sets. In Table [3] we present 16 models 
successfully derived only from sparse data that passed this test. 

Models from sparse data vs. DAM IT. Here, we used previously 
derived models based only on relative photometry which are 
stored in DAMIT. As can seen in the previous test, sparse data 
are sometimes sufficient to produce a unique model of an aster- 
oid. In 16 cases, we were able to derive a model for an asteroid 
which was already included in DAMIT and thus a model based 
on entirely different photometric data sets is available. These two 
independent models can be then compared and should be similar. 
We obtained similar resulting models for all 16 asteroids. 

Models of "mock" objects. For each asteroid shape model de- 
rived only from sparse data, we created a set of ten "mock" 
objects of roughly the same appearance and spin state (see an 
example of such shape in Fig. 0. For these synthetic objects, 
we computed their photometric data using same epochs and ge- 
ometries and with similar random noise level. These synthetic 




Fig. 2. Asteroid (810) Atossa: shape model (left panel) and an 
example of a "mock" shape model (right panel). 

photometric data sets were then used in the lightcurve inversion 
method. A check that the original model using actual data is re- 
liable is to be able to derive most of the models of the "mock" 
objects. The dispersions between the periods and pole directions 
of the "mock" objects represent the typical uncertainties of these 
parameters. For all studied asteroids, we were able to derive 
unique models for most of their "mock" objects. In all cases 
when we did not get a unique solution for the "mock" object, 
the best fit corresponded to the correct solution although other 
solutions could not be ruled out. The typical uncertainty in pole 
direction was ±10° and, for the period, ~0.1 times the difference 
between the two local minimums as determined by the period 
value and the time span of the data. 

4. An analysis of periods, poles and sizes of 
asteroid models 

Previous studies have been looking on spin state results using 
different techniques, e.g. amplitude-epoch, lightcurve inversion 
or radar methods. If there were multiple solutions for a given 
asteroid, the most probable one or simply the weighted mean 
was taken. However, this could cause systematic deviations. In 
our study, we used only results based on the lightcurve inver- 
sion method - i.e., unique and partial models presented in this 
work and models from the DAMIT. Our sample consists of 
221 asteroid models: 80 new models, 18 updates for models 
from DAMIT, 30 new partial m odels, 84 models fro m DAMIT 
and 9 new models presented bv lDurech et al.l (1201 ll) . Our sam- 
ple consists of models for 206 MB As, 10 NEAs, 3 Hungaria, 
1 Trojan, and 1 Hilda and so a statistical study is only possible 
for the MBAs. In many cases there is an ambiguity in the pole 
direction since there are two, undistinguishable mirror solutions. 
For our statistical analysis we randomly chose only one. 

In Fig. [3] we show (among other things to be discussed 
later) the relation between the proper semi-major axis and the 
proper eccentricity for asteroids in our sample and for all main 
belt asteroids. It is obvious that the positions of studied aster- 
oids strongly correlate with the MBAs population and so derived 
models are not significantly biased with respect to orbits, e.g. 
they do not lie in the inner main belt. Several asteroids in Fig. [3] 
with semi-major axis a > 3.3 AU belong to the Cybele group, 
e.g. (121) Hermione. 

4. 1 . Pole distribution analysis 

In the following study of spin axis directions, we did not use the 
Koronis family members because their spin states are correlated, 
i.e., their spin vectors are clustered towards two values of the 
obliquity dSlivanll2002l) . In Fig.|4^, we show the ecliptic latitude 
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Fig. 3. Relation between the proper semi-major axis and the proper eccentricity for asteroids in our sample and for the first 100 000 
numbered asteroids for comparison. Main resonances are shown by dotted lines. Prograde rotators are plotted with blue circles and 
retrograde rotators with red squares. The horizontal lines represent for each asteroid its estimated past drift (i.e. where the asteroid 
came from) during the collisional lifetime computed with Eq.|7] Proper elements are from the AstDyS database. 



distribution of our MBA sample. As in all similar plots, the width 
of the latitude and longitude bins corresponds to equal surfaces 
on the (A, /3)-sphere (bins are equidistant in sin j3 for latitudes 
and in A for longitudes). We confirmed the expectation that there 
is a lack of asteroids with latitudes close to the ecliptic plane. 
The latitude distribution is clearly not symmetric: about half of 
the retrograde rotators have latitudes in the bin (-53°, -90°). On 
the other hand, less than a third of asteroids with prograde spins 
are in the corresponding bin (53°, 90°). Moreover, the remaining 
prograde bins are more populated than the corresponding retro- 
grade ones. From a detailed look at the plot we can see that there 
are up to 10% more prograde rotators among the MBAs. 

In Fig. [21 we show the dependence of the ecliptic latitude /3 
of the pole direction on the diam eter D (most of the diameters 
used are based on the IRAS data dTedesco et al.ll2.QQ2h or occul- 
tations profiles and have an uncertainty +10%). Even for diam- 
eters D < 50 km, the clustering of the latitudes towards higher 
absolute values, and conversely, the depletion close to the eclip- 
tic plane is obvious and markedly so for D < 30 km. Asteroids 
with larger diameters have a more isotropic distribution of lat- 
itudes and only a moderate excess of prograde. In Fig. [6] we 
plotted the latitude and longitude distributions of asteroids with 
respect to their diameters. Based on Fig. |5] we resolved three 
different size groups: 0-30 km, 30-60 km and >60 km. 

The latitude distribution for asteroids with D > 60 km 
(Fig. [6^) is close to the uniform distribution for latitudes lower 
than 11° and for larger latitudes it exhibits an evident ex- 
cess of prograde rotators. This is in agreeme nt with theoreti- 
cal arguments presented by iDavis et alX(l 19891) and recently by 
iJohansen & Lacerdal (l2010h . On the other hand, the latitude dis- 
tribution for asteroids with D < 30 km (Fig. [(J) exhibits a strong 
depopulation of pole vectors close to the ecliptic plane (i.e. small 
absolute values of latitudes y3). The few asteroids with small lat- 
itudes have diameters D > 25 km. The latitude distribution for 
asteroids with intermediate diameters of 30-60 km (Fig. |6j;) is 
also somewhat clustered towards higher latitudes but the bins 
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Fig. 4. Fig. a) shows the ecliptic latitude distribution of all MBAs 
in our sample except the Koronis family members. The width of 
bins in latitude /3 corresponds to similar surfaces on the (A, /?)- 
sphere, so the bins are equidistant in sin/?. In Fig. b), the longi- 
tude distribution of all MBAs is plotted, again except the Koronis 
family members. The longitude bins are equidistant in A. In the 
top right corners, there are the probability values of the ^- 2 -tests 
(that the observed distributions are randomly drown from a uni- 
form distribution, see Table|2]i. 



for small latitudes are more populated. Therefor, it is probably a 
transition region between the two distinct distributions. 

It is evident that the depopulation concerns mainly objects 
with diameters D < 30 km (the distribution for the intermedi- 
ate size sample shows that the limit is probably ~50 km). This 
size roughly corresponds to the value, when the YORP effect 
starts to act and hence it is a natural candida te for a physical 
explanation. It is know n from previous studies dPravec & Harrisl 
l2000t lRubincamll2000l) that the YORP effect is significantly al- 
tering the periods and also spin vectors of these objects on a 
timescale shorter than the typical collisional lifetime of these 
objects (timescales are discussed in more details in Section 0. 
In Fig. [7] we show the relation between the spin rate and the lat- 
itude for the small (D < 30 km) and large (D > 60 km) groups 
of asteroids. In concert with YORP theory, the spin up and spin 
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Fig. 5. Dependence of the ecliptic latitude /3 (plotted as sin/?) 
of the pole direction on the models diameter D. The gray box 
shows the gap of small latitudes for asteroids with D < 30 km. 
Members of the Koronis cluster are plotted with crosses. 



Fig. 7. Dependence of the ecliptic latitude /? (plotted as sin/?) 
of the pole direction on the models spin rate, a) for asteroids 
with diameter D < 30 km and b) for asteroids with diameter 
D > 60 km. 



Observed latitude distribulion a) 
Diameter > 60 km 4.5% 



Observed longitude distribution b) 
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Observed longitude distribution d) 
30 < Diameter < 60 km 90% 
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Observed longitude distribution f) 
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Fig. 6. Histograms showing the observed latitude and longitude 
distributions of MBAs (except Koronis cluster members) for dif- 
ferent size ranges. Fig. a) shows the latitude distribution for as- 
teroids with diameters larger than 60 km, c) for asteroids with 
diameters in the range of 30-60 km, e) for asteroids with diame- 
ters smaller than 30 km, and similarly b), d) and f) for longitudes. 
The width of bins in the latitude /? and longitude lambda A cor- 
responds to similar surfaces on the (A, /?)-sphere, so the bins are 
equidistant in sin/? and A. In the top right corners, there are the 
probability values of the^-tests (that the observed distributions 
are randomly drown from a uniform distribution, see Table 23 ■ 



down and the simultaneous evolution of the latitudes towards 
higher absolute values are evident in the small asteroid sample. 

Note that observed latitude and longitude distributions can 
be biased by the convex inversion method and, therefore, this 
bias should be taken into consideration. In general, models for 
asteroids with higher amplitudes are more often successfully de- 
rived than those asteroids with lower amplitudes. This is be- 
cause, while the accuracy of the sparse photometry in both cases 
is roughly the same, the signal-to-noise ratio is significantly bet- 



ter for higher amplitude lightcurves. If we assume two bodies 
with the same shapes and orbits and different ecliptic latitudes 
of the poles, the body with higher absolute value of ecliptic lat- 
itude usually has a higher amplitude. This effect is numerically 
inv estigated in Sectionl4~3l 

ISkoglov & Eriksonl d2002h discussed the role of orbital in- 
clination in the observed distribution of latitudes. No indi- 
cations of this effect wer e found in the asteroids sample of 
iKrvszczvriska et all (120071) . We also did not find any indications 
of such correlation in our sample and so we conclude that orbital 
inclination does not affect the observed distribution of latitudes. 

We are not aware of any other physical effects in the main 
belt that could explain the non-uniform observed latitude distri- 
bution of small asteroids (D < 30 km). Collisions are believed 
to produce uniform spin distributions and close encounters with 
planets are common only among NEAs. 

There are many additional selection effects that influence the 
properties of derived models, e.g. the role of amplitude, orbit, 
time, accuracy and geometry of observations, among others. The 
significance of this bias is unknown and cannot be easily deter- 
mined. The main problem here is that for almost every asteroid, 
the photometric data are from different observers with a differ- 
ent number of measurements, quality, and purpose. The signif- 
icance of this effect can be determined only from a comparison 
of models derived from real and synthetic data of known prop- 
erties. This may be possible in a few years when the photometry 
from the Pan-STARRS is available, but not now. In the mean- 
time, the role of the selection effects seems to be small and does 
not significantly affect, for example, the latitude distribution. 

The longitude distributions of the MBAs are plotted in 
Figs. |4j) and [6}},d,f. They are, contrary to the latitude distribu- 
tions, without any statistically significant features and have very 
close to uniform distributions. The only exception are the aster- 
oids with D < 30 km, but the excess appears to be just a random 
coincidence than a result of some physical process. 

In all cases, we tested a hypothesis that the observed distribu- 
tion of latitudes or longitudes is uniform (using a^- 2 -test 9 ). The 
computed chi-squares and corresponding probabilities are listed 
in Table [2] Higher^- 2 -values and lower probabilities mean that 
the supposed hypothesis "the observed distribution is uniform" 
does not fit the observed data. If we assume a probability of 5% 
or lower as statistically significant, we can say that the latitude 
distributions for the whole sample and for asteroids with D < 30 
km and 30 < D < 60 km do not agree with a uniform distribu- 
tion. On the other hand, all longitude distributions are consistent 



9 The results based on the ^- 2 -test are also in agreement with the 
Kolmogorov-Smirnov test 
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Table 2. Test of the hypothesis that the observed pole distribu- 
tions and de-biased latitude distribution are drown from uniform 
distributions (a^f 2 -test). 
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with uniform distributions. Latitude distribution for the MBAs 
with diameters D > 60 km also disagrees with the uniform dis- 
tribution; this is because of the excess of prograde rotators. 

The overall view on the model positions within the main belt 
of asteroids, together with their estimated total drifts and the in- 
formation about whether they are prograde or retrograde rot ators 
(Fig. [3), shows behavior consistent with the Yarkovsky/YORP 
theory: there is an asymmetry of prograde and retrograde rota- 
tors near the main resonances and prograde asteroids drift out- 
wards from the Sun and can reach the resonance only from the 
left. On the other hand, retrograde rotators drift the oppos ite di- 
rection and depopulate the zone left of the resonance because 
the resonance prevents the entry of new retrograde asteroids. 
This creates an excess of prograde rotators. The same mecha- 
nism works also in the zone right of the resonance, except, in 
this case, only an excess of retrograde asteroids is now created. 
This effect is obvious in the neighborhood of the V(, and 3:1 reso- 
nances. The total drift from an asteroid's original location during 
the collisional lifetime computed with Eq. [7] are inversely pro- 
portional to the size of the asteroids. Larger asteroids (D > 50 
km) do not drift significantly while smaller asteroids frequently 
drift +0.05 AU. Note that small-sized asteroids are found mainly 
in the inner or middle part of the main belt (due to selection 
effect; they have high albedos and/or are closer to the Earth). 
Asteroids drifting through the resonances during their collisional 
lifetime are interesting from the point of the dynamical evolu- 
tion. These asteroids were either recently collisionally-affected 
or their shape models are wrong. There seems to be three such 
models, two near the 3:1 resonance and one near the resonance 
5:2. 



4.2. The Koronis family members 

The analysis of rotational state solutions for ten members of the 
Koronis aste roid family revealed a clustered distribution of their 
spin vectors dSlivarj2002tlSrivan et al.l200l. T his spi n distribu- 
tion was later explained bv lVokrouhlickv et al.l (120031) as the re- 
sult of the thermal torques and spin-orbital resonances that mod- 
ify the spin states over time of 2-3 Gyr. The modeling suggested 
an existence of two groups of asteroids: (a) low-obliquity retro- 
grade objects with rotational periods P < 5 h or P > 13 h, and (b) 
prograde rotators with periods 4 < P < 7 h that became trapped 
in a spin-orbit resonance with secular frequency S6 and thus have 
similar spin obliquities (42 -51°) and also similar ecliptic longi- 
tudes in the range of (24-73°) an d (204-259°) All t en members 
of the Koronis family studied bv lSlivanl (120021) and lSlivan et al.l 
(120031) had the expected properties: periods for pr ograde rota- 
tors w ere shifted only to higher values of 7-10 h. ISlivan et"al] 
d2009l) published spin state solutions for another four members 
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Fig. 8. Histograms of the angular differences between the gener- 
ated and derived ecliptic latitudes and longitudes of the pole di- 
rections for all successfully derived synthetic models: a) /J comp - 
ySgen, where /3 gen and p comp are generated and computed ecliptic 
latitudes and similarly b) {A comp - A gen ) cos y6 gen for ecliptic lon- 
gitudes. 



of Koronis family. Only the solution for (253) Dresda was not in 
agreement with the theoretical expectation. 

Here, we present three new models of asteroids belonging to 
the Koronis family: (832) Karin, (1482) Sebastiana, and (1635) 
Bohrmann, along with two partial models for (1350) Rosselina 
and (1389) Onnie. Only the spin state solutions for Bohrmann 
and Onnie fit the theoretical expectations. Rotational parameters 
for Karin (P = 18.3512 h, A = 242°, /3 = 46°) are outside both 
groups. Asteroids Sebastiana and Rosselina are low-obliquity 
retrograde rotators, but their periods (10.49 h for Sebastiana and 
8.14 h for Rosselina) are in the middle of the "wrong" range of 
P = 5-13 hours. Karin is the namesake and largest mem ber of 
a small and young (~ 5.8 My, iNesvornv & Botfkdl2004l) colli- 
sional family that is confined within the larger Koronis family. 
The spin state of Karin was thus likely affected during this catas- 
trophic event and changed to a random state that disagrees with 
the clustered distribution. 

We are not able to give a satisfactory explanation for 
the peculiar spin state solutions for Sebastiana and Rosselina. 
Nevertheless, we are aware of two possible scenarios: (i) the ini- 
tial rotational state and shape did not allow being captured in 
the resonance or (ii) the objects were randomly reoriented by 
non-catastrophic collisions. The timescales of such collisions 
(given by Eq. |7]i are for Sebastiana T leor ~ 7.5 Gyr and for 
Rosselina T leol ~ 14.7 Gyr. This leads to the probability of a col- 
lision during the K oronis cluster lifetime (estimated to ~2.5 Gyr, 
iBottke et al.ll200ll) ~l/3 for Sebastiana and ~l/6 for Rosselina, 
respectively, which means that random collisional reorientation 
of the spin axis is likely for at least a few of 19 asteroids in the 
Koronis cluster with known spin state solutions (most of them 
have T re or ^ 20 Gyr). 

4.3. Biases of the LI method 

We developed a numerical algorithm to estimate the selection 
effect of the lightcurve inversion method and used this approach 
to de-bias the observed distribution of asteroid's pole directions. 
The algorithm was as follows: 

1. For a model with a known shape, we randomly generated 
a new pole direction (while the overall distribution of poles 
was isotropic). 

2. For each shape with a new rotational state but with the period 
unchanged, we computed synthetic lightcurves for the same 
epochs as observed ones. 

3. To each data point i, we added the corresponding noise 
given by: 
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Fig. 9. Fig. a) shows the simulated ecliptic latitude distribution 
of all successfully derived models. In Fig. b), the de-biased ob- 
served latitude distribution is plotted. The bins are equidistant in 
sin/?. In the top right corner, there is the probability value of the 
^ 2 -test (that the observed distribution is drown from a uniform 
one, see Table[2]>. 
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Fig. 10. Histograms showing the simulated and corrected (de- 
biased) observed latitude distributions of MB As for different size 
groups, excluding the Koronis cluster members. Fig. a) shows 
the simulated latitude distribution for asteroids with the diame- 
ters larger than 60 km, c) for asteroids with the size in range of 
30-60 km, e) for asteroids with diameters smaller than 30 km. 
Similarly Figs, b), d) and f) show the observed latitude distri- 
butions corrected by the bias of the LI method. The bins are 
equidistant in sin/?. In the top right corners, there are the proba- 
bility values of the ^ 2 -tests (that the observed distributions are 
drown from a uniform distribution, see Table |2)- 



4. Finally, we performed a lightcurve inversion the same way 
as with the actual data and tried to derive a model. 

5. We repeated steps 1-4 for 50 random poles for each asteroid 
model. 

In this simulation, we used 80 models derived from com- 
bined dense and sparse data sets and 89 models from the 
DAMIT. 

For each successfully derived model we have the generated 
pole direction (A gen , /3 gen ) and period P gen , and also the com- 
puted pole and period: A comp , f3 comp , P CO mp- We computed the 
angular differences between the generated and derived ecliptic 
latitude and longitude of the pole direction: /3 comp - /3 oen and 

w 



comp 



/Igen) cos y8g en (the cos yS gen factor is used for the correc- 
tion of the different distances of meridians near the equator and 
poles). In Fig. [8] we show the histograms of these differences in 
a) ecliptic latitudes and b) ecliptic longitudes. We assumed that 
the histograms can be described by a normal distribution and we 
computed the mean and the standard deviation (fi; <f). We found 
values of (-0.2;10.2) for latitudes and (-0.2;5.2) for longitudes. 
The standard deviation <x is directly related to the typical uncer- 
tainty that we can expect in pole determination by the lightcurve 
inversion method, which is ~5°/ cos/? in A and ~10° in /?. 

In Fig. [9^, we constructed a histogram of the latitude dis- 
tribution for all successful models. The bins in j3 were again 
equally spaced in sin/?. The latitude distribution of all generated 
models was not exactly uniform, the amount of latitudes in bins 
slightly differed. To remove this effect, we divided the latitude 
distribution of successfully derived models by the latitude dis- 
tribution of all generated models normalized to unity. This cor- 
rection was also applied to latitude distributions in Figs.fTOh.c.e. 
It is obvious that the LI method is more efficient for asteroids 
with higher \/3\. The amount of successfully derived models with 
|/?| -0° is about 30% lower than with |/?| > 53°. 

In Figs. [TOh.c.e. we constructed the histograms of latitude 
distributions for successfully derived models and distinguished 
three size ranges. All three plots look very similar, except that 
with decreasing size, the ratio between models with |/J| ~0° and 
\P\ > 53° goes down. This ratio is -75% for D > 60 km, -65% 
for 30 < D < 60 km and -60% for D < 30 km. 

The histograms in Figs. [9^ and [TOk.c.e define the bias in 
latitude of the LI method and can be used for de-biasing the 
observed latitude distributions presented in Figs. |4fc and[6^,c,e. 
The de-biased histograms of latitudes are plotted in Figs.[9|3 and 
[TOb.d.f. The histograms changed only slightly and the conclu- 
sions from Section l4~T1 are still valid, i.e., the latitude distribution 
differs significantly from a uniform distribution, and especially 
so for D < 30 km. The distribution of latitudes for asteroids 
with D > 60 has an evident excess of prograde rotators while 
the distribution for a subsample with 30 km < D < 60 km shows 
an enrichment of asteroids with large latitudes (|/?| > 53°). Other 
bins have similar populations. We also performed a^ 2 -test in the 
same way as for the observed distributions (see Table [2] columns 
6 and 7). 

We did not find any significant correlation between the eclip- 
tic longitude and the efficiency of the model determination. 



Si = 



r mod 



(2) 5. A theoretical model of the latitude distribution 



where L° bs is i-th brightness observed and L™ od is ;-th bright- 
ness computed, both for the original model. This gave us 
synthetic lightcurve data equivalent to the original observed 
data, but for a new pole direction. 



In order to understand observations of main-belt asteroids, 
namely the de-biased distribution of their ecliptic latitudes [5 
(Fig. |9j)), we constructed a simple model for spin evolution 
that contains the following processes: (i) the YORP effect, 
i.e., torques arising from the emission of thermal radiation, 
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Fig. 11. Left panel: The distribution of the ecliptic latitudes B for the observed asteroids (thick line) and 100 synthetic samples 
(generated with a different random seed) at time t — (thin lines). Right panel: steady-state synthetic latitude distributions at 
time t — 4Gyr (thin lines), evolved by the YORP effect, collisions, mass shedding and spin-orbit resonances. We also applied an 
observational bias to it. The steady state was reached already within 1 Gyr. 



(ii) random reorientations induced by non-catastrophic colli- 
sions, and (iii) mass shedding after a critical rotational frequency 
is reached. 

On the other hand, we did not include gravitational torques 
of the Sun, spin-orbital resonances, damping (dissipation of ro- 
tational energy), or tumbling. Even though individual asteroids 
may be substantially affected by these processes, our model is 
for a large statistical sample of asteroids and the effect on the 
overall latitude distribution is assumed to be only minor. For 
example, gravitational torques and spin-orbital resonances usu- 
ally cause large oscillations of B for prograd e-rotating asteroids, 
but th ey remain bound to a certain interval (lVokrouhlickv et all 
l2006h . Moreover, we tried to account for these (rather random) 
oscillations in our model as well (see below). 

Our sample of 220 asteroids was the same as the observed 
sample discussed in Section [4] This means that the orbits and 
sizes correspond to re al asteroids . The m odel for spin evolution 
was similar to that of iBroz et al. where it was used for 

studies of the long-term evolution of asteroid families. We as- 
sumed the following relations for the rate of the angular veloc- 
ity to and the obliquity e due to the YORP effect 



Tt =fM ' 



i = 1 ...200. 



At 



gi(e) 



(3) 
(4) 



where /- and g-f unctions were given bv lCapek & Vokrouhlickvl 
d2004l) for a set of 200 shapes with mean radius Rq — 1 km, bulk 
density po = 2500 kg/m 3 , located on a circular orbit with semi- 
major axis «o = 2.5 AU. We assigned one of the artificial shapes 
(denoted by the index ;) randomly to each individual asteroid. 10 
We had only to scale the /- and ^-functions by a factor 



cyorp 



ao 



-2 i v-1 
Pbulk I 

P0 / 



(5) 



10 We did not use the convex-hull shapes derived in this work for two 
reasons: (i) the two samples of shapes are believed to be statistically 
equivalent and it is thus not necessary to compute the YORP torques 
agai n; (ii) the YORP effect seems sensitive to small-scale surface struc- 
ture llScheeres & Mirr ahimi 2007) which cannot be caught by our shape 
model. Nevertheless, the YORP torque remains of the same order, so the 
random assignment of shapes seems reasonable. 



where a, R, pbulk are semi-major axis, radius, and density of the 
simulated body, respectively, and cyorp is a free scaling parame- 
ter, which can account for an additional uncertainty of the YORP 
model. Because the values of f's and g's were computed for 
only a limited set of obliquities (with a step Ae _ 3 0°) we used 
interpolation by Hermite polynomials (lHilllll982l) of the data 
in ICapek & Vokrouhlickvl (120041) to obtain a smooth analytical 
functions for f )(e) andg,(e). 

When the angular velocity approached a critical value (i.e., 
the gravity was equal to the centrifugal force) 



(6) 



we assumed a mass shedding event. We kept the orientation 
of the spin axis and the sense of rotation but reset the or- 
bital period P = 2tt/lo to a random value from the interval 
(P\,Pi) = (2.5,9) hours. We also altered the assigned shape 
since any change of shape can produce a different YORP effect. 
We did not change the mass, however. 

The differential equations (01 and (|4]i were integrated nu- 
merically by a simple Euler integrator. The usual time step was 
At = lOOOyr. The time scale of the spin axis evolution for small 
bodies (D 10 km) is tyorp - 500 Myr. After 3 times tyorp 
most of these bodies have spin axes perpendicular to the ecliptic. 

We also included a Monte-Carlo model for spin axis reorien- 
tations caused by collisions. 11 We used an estimate of the time 
scale bv lFarinella et al.1 d!998l) 



B [^) (£) 



(7) 



where B = 84.5 kyr, B\ = 5/6, B 2 = 4/3, D = 2m and w cor- 
responds to period P = 5 hours. These values are characteristic 
for the main belt. After a collision, we reset the spin axis periods 
to random values, using the interval (P' v P' 2 ) = (2.5, 9) hours for 
the period. Since the time scale is T reol 3 Gyr for the small- 
est (D 5 km) bodies, reorientations are only of minor impor- 
tance. However, note that the probability of the reorientation is 
enhanced when the YORP effect drives the angular velocity u> 
close to zero. 



11 Collisional disruptions are not important, since we are only inter- 
ested in the steady state. We can imagine that whenever an asteroid from 
our sample is disrupted, another one with a randomly oriented spin axis 
is created by a disruption of a larger body. 
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There were several free parameters in our model: the cyorp 
parameter, thermal conductivity K, bulk density pbuik, initial dis- 
tribution of (3 and initial distribution of a>. 

Our aim was to start with a simple /?- and w-distribution, wait 
until a steady state was reached, and then compare the resulting 
synthetic to observed latitude distributions. We applied an obser- 
vational bias derived in Section l4~3l to the synthetic distribution. 

We partly accounted for spin-orbital resonances acting on 
prograde asteroids by adding a sinusoidal oscillations to /? with 
a random phase and an amplitude 40°, which are typically 
induced by resonances. This procedure naturally decreased the 
right-most bin (sin/? = (0.8, 1)) of the synthetic distribution and 
increased the next bin (sin/? = (0.6, 0.8)). 

We started with reasonable parameters of cyorp = 0.33, K = 
10~ 2 W/K/m, pbuik = 2500 kg/m 3 , a Maxwellian distribution 
of to, a uniform distribution of sin f3 (i.e., an isotropic distribution 
of spin axes). We ran 100 such simulations with different random 
seeds. A steady state was reached within 1 Gyr. The resulting 
latitude distributions are shown in Figures QT|and Q~2] 

From these it can be seen that: (i) the observed distribution 
of /3 for small asteroids seems compatible with our model; the 
YORP effect is capable of creating such an uneven distribution 
and (ii) there is a discrepancy for large asteroids (especially in 
bins sin/? e (-1, -0.8) and (0.2, 0.4)), which can be explained as 
a pre ference for prograde rotators in the primordia l population 
(see iDavis et al.ll 19891 : iJohansen & Lacerdall2010l). T he results 
regarding the spin rates agree with lPravec et al. I (12008). so we do 
not repeat the discussion here. 

We also tested the sensitivity of our results with respect 
to the free parameters. The thermal conductivity did not seem 
important (we tested K = 10~ 3 W/K/m). A simulation with 
cyorp = 0.66, p = 1300 kg/m 3 , and a uniform distribution of 
orbital periods P oc 1/cj produced almost the same resulting lat- 
itude distribution. Nevertheless, a value of cyorp = 1 -00 seems 
too high because the extreme bins of the /^-distribution were 
overpopulated. 12 There is only a weak dependence of our results 
on the period ranges that we used for resetting the orbital period 
after a mass-shedding event (P\,P2) and collision (P' V P' 2 ). As 
would be expected, the values P\,P2, P[,P'o significantly affect 
the period distribution. The relatively weak dependence on the 
free parameters likely stems from the fact that we presume the 
steady-state. Even though the free parameters change, e.g., the 
strength of the YORP effect and the evolution of spins is slower 
and/or faster, after reaching a steady state, the basic characteris- 
tics of the latitude distribution remain similar. The observed /?- 
distribution of small asteroids (D < 30 km) cannot be explained 
by our simulation without accounting the YORP effect. 

6. Conclusions 

The results of this paper can be summarized as follows. 

We used combined dense and sparse data to derive new aster- 
oid shape models. We systematically gathered and processed all 
available sparse photometry from astrometric surveys and em- 
ployed valuable data from seven observatories (see Table [1} in 
lightcurve inversion. 

We derived 80 new unique models of asteroids, from which 
16 are based only on sparse data. We also present 30 partial mod- 
els with accurate rotational periods and estimated ecliptic lati- 

12 As an alternative hypothesis, we assumed the spin axis evolu- 
tion without a YORP effect (cyorp = 0). In this case, the initial /?- 
distributions (Fig. fTT] left panel, thin lines) do not change significantly 
in time. 



D > 60 km 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 
sin p 



Fig. 12. Steady-state synthetic latitude distributions for three dif- 
ferent size ranges: D > 60 km (top), 30 km < D < 60 km (mid- 
dle) and D < 30 km (bottom). The synthetic distributions are 
plotted by thin gray lines while the observed distribution by a 
thick black line. 

tudes of the pole directions and 18 updated solutions based on 
new data for asteroids already included in DAMIT. 

In the future, quality sparse data sets will be produced by all- 
sky surveys such as Pan-STARRS, the Large Synoptic Survey 
Telescope (LSST), and the Gaia satellite. When these data are 
available, we will be able to apply the same methods to in order 
to derive many more new asteroid models. These surveys will 
have one advantage over dense data: their selection effects (e.g., 
with respect to the orbit) will be known. This will allow us to 
make a more accurate analysis of the asteroid population. 

As expected, the observed ecliptic longitude distribution of 
asteroid spin vector is independent of diameter and is compati- 
ble with a uniform distribution. Unlike the latitude distribution, 
the observed ecliptic longitude distribution is not significantly 
biased by the LI method. However, the effect of the LI bias is 
only minor and the global features of the observed latitude dis- 
tribution do not change. The observed (and de-biased) latitude 
distribution for asteroids with diameters D > 60 km shows an 
excess of prograde rotators in the latitude interval (11°, 90°). 
This excess is probably primordial. On the other hand, the lati- 
tude distributions for the entire sample and in particular for as- 
teroids with D < 30 km, is strongly anisotropic. 

The dynamical evolution of asteroid spins seems to be dom- 
inated by the YORP effect and also by collisions and mass shed- 
ding for asteroids with diameters D < 30 km. We calculate that 
YORP (with a small contribution for the LI method's bias) is ca- 
pable of producing the observed depopulation of spin vectors for 
small asteroids. 

We are not yet able to study small asteroids in individual 
families (small bodies at the outskirts of a family should have 
extreme spins); this is an aim of future work. 
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Table 3. List of new asteroid models derived from combined data sets or sparse data alone. 
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Table 3. continued. 



J. Hanus et al.: A study of asteroid pole-latitude distribution 



Asteroid A { X 2 j8 2 P N lc N m N 6m N 699 N 103 N El2 N m6 N 950 N Wp 







[deg] 


[deg] 


[deg] 


[deg] 


[hours] 












1482 


Sebastiana 


262 


-68 


91 


-67 


10.48965 


2 


1 131 


39 


30 




1514 


Ricouxa 


251 


75 


68 


69 


10.42468 


3 


1 68 


56 






1568 


Aisleen 


109 


-68 






6.67597 




82 


37 






1635 


B ohrmann 


5 


-38 


185 


-36 


5.86427 


8 


1 108 


47 






1659 


Punkaharju 


259 


-71 


75 


-22 


5.01327 


2 


1 118 


66 






1682 


Karel 


232 


32 


51 


41 


3 37485 




54 


84 




36 


1709 


Ukraina 


165 


-61 


2 


-40 


7.30517 


2 


1 46 


79 






1742 


Schaifers 


198 


57 


47 


55 


8.53270 


3 


1 106 








1747 


Wright 


227 


31 






5.28796 




70 


55 






1889 


Pakhmutova 


22 


-76 


167 


-40 


17.5157 




68 


46 


35 




1930 


Lucifer 


32 


17 


211 


-19 


13.0536 


6 


1 106 


43 


66 




2156 


Kate 


49 


74 






5.62215 


4 




44 






3678 


Mongmanwai 


125 


-65 






4.18297 


2 




103 




31 


4483 


Petofi 


107 


40 






4.33299 


3 




36 







Notes. For each asteroid, the table gives also the number of dense lightcurves N lc observed during N app apparitions and the number of sparse data 
points for the corresponding observatory: N^gg, A^99, Neu, Ng96, Miso and A^Hi p . Pole solutions preferred by asteroid occultation measurements 
(Dur ech et al.l201lh are emphasized by a bold font. 
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J. Harms et al.: A study of asteroid pole-latitude distribution 
Table 4. List of partial models derived from combined data sets. 



Asteroid A P N Ic N app N 689 N 699 N 703 N E12 N G96 N 950 N Hip 







[Cleg] 


Meg] 


[hours] 
















163 


Erigone 


-60 


14 


16.1403 


3 


1 


168 




72 






187 


Lamberta 


-58 


9 


10.66703 


9 


2 


159 




52 


53 




233 


Asterope 


49 


8 


19.6981 


13 


3 


184 




80 




165 


272 


Antonia 


-70 


6 


3.85480 


5 


1 


109 




60 




36 


281 


Lucre tia 


-54 


11 


4.349710 


6 


3 


123 


30 


62 






313 


Chaldaea 


33 


18 


8.38992 


9 


3 


176 




80 






390 


Alma 


-60 


22 


3.74116 


2 


1 


109 






34 




510 


Mabella 


-59 


12 


19.4304 


6 


2 


145 




60 






550 


Senta 


-63 


13 


20.5726 


9 


1 


151 






61 




622 


Esther 


-61 


9 


47.5042 


5 


1 


120 




60 






692 


Hippodamia 


-52 


25 


8.99690 


3 


1 


1 14 




78 


32 




733 


Mocia 


36 


16 


11.37611 


2 


1 


175 




44 






746 


Marlu 


-54 


18 


7.78887 


3 


1 


133 




47 


34 




784 


Pickeringia 


58 


15 


13.1699 


1 


1 


188 




67 


32 




823 


Sisigambis 


57 


9 


146.58 


8 


1 


123 




90 






877 


Walkure 


53 


12 


17.4217 


3 


1 


141 


45 


104 




32 


899 


Jokaste 


-58 


19 


6.24812 


3 


1 


140 






43 




1010 


Marlene 


46 


7 


31.066 


8 


1 


104 




52 






1103 


Sequoia 


-48 


19 


3.037977 


2 


1 


111 




36 


30 




1185 


Nikko 


46 


12 


3.78615 


3 


1 


91 




46 


32 




1188 


Gothlandia 


-63 


19 


3.491820 


2 


1 


129 


33 


67 


41 




1214 


Richilde 


-59 


15 


9.86687 


4 


1 


101 




78 






1282 


Utopia 


-39 


21 


13.6228 


4 


1 


116 




72 






1350 


Rosselia 


-58 


13 


8.14011 


1 


1 


114 




48 






1368 


Numidia 


-50 


14 


3.640740 


3 


1 


129 




47 






1379 


Lomonosowa 


-62 


17 


24.4845 


2 


1 


96 




100 






1389 


Onnie 


-56 


10 


23.0447 


2 


1 


85 


33 


47 


32 


40 


1665 


Gaby 


49 


17 


67.905 


1 


1 


81 




80 






1719 


Jens 


-56 


19 


5.87016 


2 


1 


78 




48 


40 




2001 


Einstein 


-51 


22 


5.4850 


2 


1 






84 







Notes. For each asteroid, there is the mean ecliptic latitude f3 of the pole direction and its dispersion A, the other parameters have the same meaning 
as in Table [3] 



Table 5. List of improved asteroid models that were originally designated in DAMIT as "preliminary". 



Asteroid 






Ai 


Pi 


P 


Ni c 


N 

1 ~ app 


^689 


N m 


N 10 3 


N E u 


N G 96 


A 95 o AW 




[deg] 


[deg] 


[deg] 


[deg] 


[hours] 


















73 Klytia 


266 


68 


44 


83 


8.28307 


21 


7 


131 


36 


98 


47 






82 Alkmene 


164 


-28 


349 


-33 


13.00079 


11 


1 


158 




72 


36 


38 


192 


132 Aethra 


326 


67 






5.16827 


4 


2 


204 




55 








152 Atala 


347 


46 


199 


61 


6.24472 


2 


1 


101 




32 








277 Elvira 


121 


-84 






29.6922 


22 


5 


142 




36 


51 






278 Paulina 


307 


31 


118 


38 


6.49387 


3 


1 


195 




51 








311 Claudia 


214 


43 


30 


40 


7.53138 


23 


6 


114 


33 


108 




40 




484 Pittsburghia 


70 


46 






10.64977 


2 


1 


100 




52 








516 Amherstia 


254 


22 


81 


54 


7.48431 


5 


3 


162 




32 








534 Nassovia 


66 


41 


252 


42 


9.46889 


16 


6 


151 




64 


32 






614 Pia 


348 


48 


162 


27 


4.57872 


2 


1 


121 




78 








714 Ulula 


224 


-10 


41 


-5 


6.99837 


9 


2 


177 




67 








770 Bali 


70 


50 


262 


45 


5.81894 


2 


1 


131 




52 








915 Cosette 


350 


56 


189 


61 


4.469742 


1 


1 


106 




32 


35 






1012 Sarema 


45 


67 


253 


63 


10.30708 


2 


1 


74 




42 








1022 Olympiada 


46 


10 


242 


52 


3.83359 


5 


2 


107 




91 








1088 Mitaka 


280 


-71 






3.035378 


1 


1 


104 




39 


41 






1223 Neckar 


252 


28 


69 


30 


7.82401 


16 


7 


132 


33 


60 




36 





Notes. Pole solutions preferred by asteroid occultation measurements (Durech et al. 201 111) are emphasized by a bold font. 
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J. Harms et al.: A study of asteroid pole-latitude distribution 



Table 6. Observations used for successful model determination that are not included in the UAPC. 



Asteroid 


Date 




(lh c p rvpr 




A etpmi c\ 
rvs lci vjiu 


Date 




(ih q p rvpr 


13 


Egeria 


2007 


/ \ — * a a a o 

9-2009.3 


p% " 1 1 /^AAA\ 

Prlcher (2009) 


808 


Merxia 


2003 


1 26.8 


Casulli" 


14 


Irene 


2007 


1 i— 2009.5 


Pilcher (2009 J 






2003 


1 28.8 


Casullr 






2008 


2 3.7 


Polishook 






2003 


2 7.0 


Bernasconi" 


40 


Harmonia 


2008 


1 ^ / ~W \ A A 1 

12-2009.1 


T"l"1 1 AAA \ 

Pilcher (2009) 






2003 


2 8.0 


Bernasconi" 






2010 


5-2010.6 


Prlcher (2010) 


832 


Karin 


2003 


O TAAn A 

8-2003.9 


Yoshida et al. (2004) 


68 


Leto 


2008 


1 26.0 


Pilcher* 






2004 


f\ A A A A 

9-2004.9 


Ito & Yoshida (2007) 






2008 


1 30.0 


Pilcher* 


899 


Jokaste 


2003 


11-2003.12 


Stephens (2004) 






2008 


2 6.1 


Pilcher* 


1010 


Marlene 


2005 


1-2005.3 


Warner (2005 a) 


264 


Libussa 


2005 


2-2005.3 


Pilcher & Coonev (2006) 


1022 


Olympiada 


1999 


6-1999.6 


Warner (2005b) 






2008 


10-2008.12 


Pilcher & Jardine (2009) 


1087 


Arabis 


2003 


2 23.0 


Lehky 


272 


Antonia 


2007 


12-2008.1 


Pilcher (2008) 


1140 


Crimea 


2005 


4-2005.4 


Stephens (2005a) 


310 


Margarita 


2010 


3-2010.5 


Pilcher & Oey (2010) 


1185 


Nikko 


2004 


11-2004.11 


Stephens (2005b) 


390 


Alma 


2004 


8-2004.8 


Stephens (2005c) 


1282 


Utopia 


2000 


11-2000.11 


Warner (2011b) 


400 


Ducrosa 


2005 


1-2005.1 


Warner (2005a) 


1333 


Cevenola 


2002 


2-2002.2 


Warner (2002) 


436 


Patricia 


2U02 


12—2003.1 


Wamer ( 2UU3 ) 


1635 


Bohrmann 


2003 


9—2003.10 


Stephens & Wainei (2UU4) 


544 


Jetta 


2004 


8-2004.8 


Stephens (2005c) 


1659 


Punkaharju 


2000 


11-2000.11 


Warner (2011c) 


573 


Recha 


2001 


1-2001.1 


Warner (201 lb) 


1719 


Jens 


2000 


9-2000.9 


Warner (2011c) 


714 


Ulula 


2005 


9 23.8 


Henych 


1930 


Lucifer 


2003 


10-2003.10 


Warner (2005a) 






2005 


9 25.8 


Henych 


2001 


Einstein 


2004 


12-2004.12 


Warner (2005a) 






2005 


9 30.8 


Henych 


3678 


Mongmanwai 


2003 


3-2003.3 


Stephens (2003) 






2005 
2005 


10 1.8 
10 10.8 


Henych 
Henych 











Notes. (fl) On line at http://obswww.unige .ch/~behrend/page_cou.html ( ' On line at http://aslc-rvm.org/Pilcher. html 
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